RUB-TPII-2/94 
May, 94 



Electromagnetic Nucleon Properties and 
Quark Sea Polarization 

Chr.V. Christov I'SiA.Z. GorskiBlK. GoekeB§and P.V. Pobylitsai 

^ Institut fiir Theoretische Physik II, Ruhr-Universitdt Bochum, D-44780 
Bochum, Germany 

^ Institute of Nuclear Physics, Radzikowskiego 152, 31-342 Cracow, Poland 

^ Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg 188350, 
Russia 



Abstract: In this paper we present the derivation as well as the numerical re- 
sults for all electromagnetic form factors of the nucleon within the chiral quark 
soliton model in the semiclassical quantization scheme. The model is based on 
a semibosonized Nambu - Jona-Lasinio lagrangean where the boson fields are 
treated as classical ones. Other observables, namely the nucleon mean squared 
radii, the magnetic moments and the nucleon-A splitting are calculated as 
well. The calculations have been done taking into account the quark sea po- 
larization effects. The final results, including rotational 1/Nc corrections, are 
compared with the existing experimental data and they are found to be in a 
good agreement for the constituent quark mass of 400 — 420 MeV. 



Permanent addressdnstitute for Nuclear Research and Nuclear Energy, Sofia, 
Bulgaria 

^ christov@neutron.tp2.ruhr-uni-bochum.de 

^ gorski@bron.ifj.edu.pl 

^ goeke@hadron.tp2.ruhr-uni-bochum.de 



To appear in Physics Reports (1995) 



1 Introduction 



In the last years one of the most challenging problems in elementary particle 
physics seems to be the solution of QCD in the low energy region. The main 
difficulties are due to the non-perturbative effects caused by the growing ef- 
fective coupling constant of the fundamental theory in the low energy limit. 
This prevents one from using the well-known main tool of theoretical physics 
— the perturbation theory. Because of this, the most intriguing features of 
QCD - confinement and chiral symmetry breaking - still remain conceptual 
and practical problems. 

The above mentioned obstacles have initiated an increasing interest among 
the physicists in non-perturbative methods and effective low-energy models of 
hadrons. The effective models are expected to mimic the behaviour of QCD 
at energies below ~ 1 GeV (confinement and/or chiral symmetry breaking) 
and to reproduce experimental data in this region. In principle, these models 
could be related to QCD by integrating out gluonic fields and reparameterizing 
fermionic degrees of freedom. 

The simplest purely fermionic Lorentz invariant model with spontaneous chiral 
symmetry breaking is the Nambu-Jona-Lasinio (NJL) model j^. It contains 
chirally invariant local four-fermion interaction terms. The NJL Lagrangian 
in its simplest SU(2) form has the following structure: 

Cmjl = ^{t^-mo)^ + ^ [{^^y + {mr^.^y ] , (1) 

where G is the coupling constant, mo is the current quark mass and r are the 
Pauli matrices in the isospin space. We consider the up and down quarks to 
be degenerated in mass. 

The NJL model is generally solved after applying the well known bosonization 
procedure following Eguchi to arrive at the model expressed in terms of 
auxiliary meson fields a, vr: 

Cnjl = ^ [i^-mo- g^{a + zvf ■ T^r^) ] ^ - ^{a^ + vf^). (2) 

The new couphng constant /i^ is related to the initial G hj G = gl/fJ^"^- Here, 
(/tt is the physical pion-quark coupling constant implying that vr is the physical 
pion field. The meson fields are constrained to the chiral circle: 

a\x) + n\x) = , (3) 

where = 93 MeV is the pion decay constant. 

In the chiral quark soliton model l^-Q based on the lagrangean (frequently 
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referred simply as NJL model) the baryons appear as a bound state of Nc va- 
lence quarks coupled to the polarized Dirac sea. Operationally the baryon 
sector of the model is solved in two steps. In the first step, motivated by 
the large limit, a static localized solution (soliton) is found by solving 
the corresponding equations of motion in an iterative self-consistent pro- 
cedure assuming that the a and vf fields have hedgehog structure. However, 
this hedgehog soliton does not preserve the spin and isospin. In order to de- 
scribe the nucleon properties one needs nucleon states with the proper spin 
and isospin numbers. To this end making use of the rotational zero modes 
the soliton is quantized. In fact, such a cranking procedure was elaborated in 
refs. and problems arising in this scheme due to the regularization are 
considered at first in ref. [0. Within the semiclassical quantization scheme 
quite successful calculations including Dirac sea polarization effects have been 
reported for the nucleon-delta mass splitting magnetic moments and 
axial-vector form factor [|1^ as well as some results for the nucleon electric 



form factors |TT|. Very recently, in the semiclassical quantization procedure 



important l/N^ rotational corrections have been derived |T2| which improve 



considerably the theoretical values for the isovector magnetic moment. 

Because of several reasons the model, described by lagrangean (H), is con- 
sidered as one of the most promising effective theories describing low energy 
QCD phenomena. 

First, the model is the simplest quark model which provides mechanism for 
spontaneous breaking of the chiral symmetry — the basic feature of QCD. The 
mesons appear as composite states of quark- antiquark pairs. The philosophy 
behind this approach is based on the hypothesis that the chiral symmetry 
breaking is the dominant mechanism which determines the structure of the 
low-mass hadrons whereas the confinement has no direct impact om them |Q] . 

Second, the self-consistent solitonic solution has been found to exist in the 
large Nc limit of the model within the physically acceptable range of parame- 
ters fixed in the meson sector. The chiral soliton provides a good description 
of the nucleon and gives us possibility to take into account vacuum polariza- 
tion effects from the Dirac quark sea. Calculations done in the semiclassical 
quantization procedure yield quite reasonable results for quantities like the 
nucleon radii, S-terms, axial vector coupling constant, A-nucleon splitting as 
well as the mass splittings within the octet and decuplet. 



Third, there are various hints [|l^-[T5| that the NJL-type Lagrangian can be ob- 
tained from QCD in various low-energy approximations. It should be stressed 
that the large Nc limit plays a prominent role in those considerations. 

The aim of this paper is to calculate the electromagnetic nucleon form factors 
within the SU{2) chiral quark soliton model based on the semibosonized NJL- 
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type lagrangean with the vacuum polarization effects taken into account. The 
calculations done in the semiclassical quantization scheme will include the 
rotational l/N^ corrections, which have been shown |T^ to be important for 



the isovector magnetic moment. Calculations of this sort give us possibility to 
calculate the nucleon electromagnetic form factors, as well as such quantities 
like the electric mean squared radii, the magnetic moments, the nucleon-A 
energy splitting, and the electric and magnetic charge distributions of the 
nucleon. 

The paper is organized as follows: we begin in Section 2 with the electromag- 
netic current in the model defined in terms of path integrals and introduce the 
rotational zero modes treating the angular velocity as a perturbation. After 
that we consider the problem of regularization. In Section 3 we derive the ex- 
pressions for the form factors including terms up to the linear order in angular 
velocity {1/Nc). In Section 4 we present and discuss our numerical results. 



2 The current and regularized action 



Our main goal in this section is to introduce the matrix element of the elec- 
tromagnetic currents in terms of path integrals and to evaluate it within the 
semiclassical quantization scheme. We follow the valence picture for the nu- 
cleon in which it appears as a bound state of N^. valence quarks coupled to 
the Dirac sea. Since the NJL model is not renormalizable a special attention 
will be paid to the problem of regularization. 

We start with the definitions of the electromagnetic current of a fermion field 

J^{x) = ^^(x)7o7MQ^(a^) • (4) 
Here Q is the quark charge matrix 

Q = ii + lr3. (5) 
^62 ^ ' 

Using the partition function of the model in Minkowski space 

we express the nucleon matrix element of current as a path integral: 

{N{p/) I J^(0) I N{p)) = lim - / d^x d^ e'P''^ e'^P'^ f VU f D^P^^ 
xJ;v(T/2,f)4(-r/2,f)v[/t(o)7o7;.g^(0) e^U'^^'^(^)^ . (J) 
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written in terms of quark \I/ , and meson U fields. The equality (J^) should 
be understood as a limit at large Euclidean time separation. Here Z is the 
normalization factor which is related to the same path integral but without 
the electromagnetic current J^. In fact, the latter represents the correlation 
function of two nucleon currents Jn{T/2,x). D is the operator 

D{U)=tdt-h{U), (8) 

which includes the one-particle hamiltonian 

a ■ V 

h(U) = + PMU^' + moP . (9) 

^ 

Here a and (5 are the Dirac matrices and the constituent quark mass M = gf^^. 
The meson fields can be equivalently written as 

U= e^-^P^^ (10) 

where 

— * 

TT = — and |7f| = /7rSin(P(r)) (11) 

Kl 

It can be easily checked that the hamiltonian is hermitian: /i^ = h. The com- 
posite Nc quark operator with the quantum numbers JJ^^TT^ (spin, isospin) 
of the nucleon can be chosen to be one of the loffe currents ilH 

Jn{x) = ^e^--''-^ri53;/-;>vl/^^^^(a;) . ■■^p,j,^{x), (12) 

where is the color index, and Fjj^ is a matrix in /j standing for both 
the flavor and the spin indices. 

In eq.(^ the quarks can be integrated out: 

{NmUO)\N{p)) = ^ N^ryi^Tj^jtrn J d'xd^y e^'- e-^^^ J VV 

X I (r/2, f I I 0, 0) f,f, (7o7mQ) fg' (0, | | - T/2, y}^,,, - Tr ((0, | | 0, 0) 7o7mQ) 

x(T/2,x|-| -T/2,y)f,,,jll{T/2,x\-\ -T/2,y)f^,^ e^'^ff , (13) 

where the effective action 

Seff = -zN^Tilog D{U). (14) 

Since we treat the meson fields classically {i.e. at the 0-loop level) and also 
they are constrained on the chiral circle (0), the only non-trivial part of the 
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action is the fermionic part Trlog D{U). In eq.(|T3D in a natural way the cur- 
rent matrix element is split in a valence - the first term, and a Dirac sea 
contribution - the second one (see the diagrams on the l.h.s. of Fig.|l] a) and 
b)). Since the model is not renormalizable, in principle, the contribution of 
the Dirac sea includes divergences and should be regularized. The second part 
of this section will be devoted to the problem of regularization. 

In eq. (0) we are left with the integral over the meson fields U. In leading order 
in Nc it can be done in a saddle point approximation. To that end we look for 
a stationary localized meson configuration (soliton) of hedgehog structure 

U{x) = , (15) 

which minimizes the effective action. As it was mentioned before, the hedgehog 
soliton U {x) does not preserve the spin and isospin. As a next step, making 
use of the rotational zero modes we quantize it. Since the fluctuations which 
correspond to the zero modes are not small, they have to be treated "exactly" 
in the meaning of path integral. It can be done assuming a rotating meson 
hedgehog fields of the form 

U{x,t) = R{t)U{x)R^{t) , (16) 

with R{t) being a time- dependent rotation SU(2) matrix in the isospin space 
and 

R^R = i. 

It is easy to see that for such an Ansatz one can transform the effective action 

TrlogD(f/) = Trlog(D(f/) - 1]) (17) 
in order to separate the angular velocity matrix: 

^ = -iR+{t)R{t) = Ulara, (18) 

from D{U). Similar to the effective action the quark propagator in the back- 
ground meson fields U can be rewritten as 

{T/2,x\^^\-T/2,y) = R(T/2) {T/2, x\ ^^Jy—^ \ -T/2,y) R+i-T/2) . 

(19) 

Here, the operator D{U)—Q corresponds to the body-fixed frame of the soliton 
in which the quark fields are transformed as 

^ — > i?(t)^ and — > ^^i?^(t). (20) 

Since (as can be seen below from the canonical quantization rules) Q ^ ^ one 



6 



can consider Q as perturbation and evaluate any observable as a perturbation 
series in Q which is actually an expansion in 

In this scheme the matrix element (0) of the current can be written as 

{N{p/)\ J^{0)\N{p)) = J d^xd^ydh e'^'"^ Q-'^'^ e-'^-^'^ 

,JvRD%j,^[R{T/2)]D%j^^[R{-T/2)]{T/2,x\^^ \ -T/2,y-)^-^ 

X I (T/2, X I I 0, i) R+{Ohol,QR{0) (0, z\ I - T/2, y) 

- Tr((0, ^\^^ I 0, i?+(0)7o7MQ^(0)) (T/2, f | | - T/2, ^-)| 

Here, the finite rotation matrix (Wigner D-function) D^^fj^^j^, which carries 
the spin and isospin quantum numbers of the nucleon, appears due to the 
rotations R{t) of the valence quark propagators in eq.(|l9D correlated by the 
FjjI^ tts matrices. The D-functions represent the collective part of the nucleon 
wave function in the semiclassical quantization scheme. The integral over z is 
due to the translational zero modes treated in the leading order. 

Now we are ready to make an expansion in Q. For the effective action it 
yields up to the second order in fl: 

S,ff ^ iV.TrlogD + dtnl (22) 

Here B is the moment of inertia. However, this quantity includes a divergent 
contribution from the Dirac sea which should be regularized. Later we will 
come back to this point. A derivation of the moment of inertia including 
regularization can be found in ref. and it has been calculated numerically in 
ref. P,^. The first term in eq.(p2D will be absorbed in Z whereas the second one 
drives the evolution in the space of matrix R. Expanding the quark propagator 

^ 1 + 1^1 + ... (23) 



D-Q D D D 

we can separate the zero order (~ A*"^) and the linear order (~ ]^) corrections 
in VL. The expansion in VL is illustrated in Fig.|I] a) and b) for the valence 
contribution and for the Dirac sea one, respectively. 

After the expansion in VL we should deal with the path integral over R. In the 
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case of the zero order term we have 

R2=R{T/2) 

J dRidR2DiT3j3(Ri)Di*T3j3(R2) / VR (r+{0)QR{0))^^ 

Ri=R{-T/2) 

^ gi|/dtU^(t')^ (24) 

whereas for the hnear in Q terms we are left with a more comphcated integrah 

R2=R(T/2) 

J dRidR2DiT3j3(Ri)DiV3j3(R2) / VR [R+{0)QR{0))^^Q,{t) 

Ri=R(-T/2) 

^ g4/dt'f7,^(t')_ (25) 

For the isoscalar part |l of the charge matrix Q the integrals (0) and ( ]25|) 
reduce simply to the normalization of the finite rotation matrix Z^i'^^j^. In the 

3 

case of the isovector part ^ we use the identity 

{r^O)t'^R{0))^^ = I TT{R^O)T-R{Oy) {r')f, (26) 

in order to separate the R{t) dependent part of the current which does not 
carry flavor indices f,g. The path integrals (p^ and ( p5| ) can be taken rigor- 
ously [^. We essentially made use of the basic feature of the path integral: 

<?2=g(T2) 

j Vq FMh)) ■ ■ ■ Fn{q{tn)) e'' = (g2, ^ \ T{A(g(ti)) ■ ■ ■ A(g(tn)} I gi, Ti), 

<ii=giTi) 

namely that the path integral, which contains Fi{q{ti)), can be equivalently 
written as a matrix element of the time ordered product T of the corresponding 
operators Fi{q(ti)) . In our case we have the well-known canonical quantization 
rule 

^ § , (27) 

where Ja is the spin operator. For the zero order the collective path integral 
(p^ ) is reduced to an ordinary integral of three Wigner D-functions whereas 
the final result for (pSj) is a time ordered product: 

^{-t)Dab{Rmjc + mJcDabiRm, (28) 

sandwiched between the nucleon rotational wave functions. Now using the 
standard spectral representation of the quark propagator 

^ en>0 
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+i^{t - t') e-""(*'-*) ) $t (x) , (29) 

en<0 



it is straightforward to evaluate the matrix element of the current - eq. (^Tj) . 
Here and e„ are the eigenfunctions and eigenvalues of the single-particle 
hamiltonian h: 

+ (3MW^ + mo/3)$„, = e„$„, . (30) 

Now we come to the problem of regularization. Since the NJL model is a non- 
renormalizable theory, a regularization scheme of an appropriate cut-off A is 
needed to make the theory finite. In our case it concerns the divergencies in the 
Dirac sea contribution to the matrix element of the electromagnetic current 
- diagrams in Fig.|I]b). In Minkowski space-time it can be equivalently written 
as: 

{N{pi)\ J^(0)|iV(p-))_ = zi j d^z e-(p'-P>^Y^^^D-T3J3[R(T/2)] 

(31) 
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Here in the effective action S^ff (0) we include an explicit electromagnetic 
coupling in a minimal way: 

D{U) = idt - h{U) - n- Q^ol^.A^ , (32) 

where is an external electromagnetic field. 

We will work in Euclidean space-time where we use the metric tensor's sig- 
nature transformed from (H ) to ( ), i.e. the Euclidean metric 

tensor g^^ = — 5^,^ and the general formulae to perform the transformation 
from Minkowski to Euclidean space are: 

A* _ „• /lO Ak _ Ak 

~ ''^M 1 ^E ~ 1 

lE = ilM 1e = 1m. 

T = X E "^-^M ~ ' "^E "^M ' 

J d^XE = i J d^XM , dr = -idt , 

{7?,7f} = -2V, {7f,75^} = 0, 

(7?)^ = -7? , (75^)^ = 75^ = 7^ ■ (33) 



The transformation of the electromagnetic field under the Wick rotation 
(transformation to the Euclidean space-time) is defined in (^) and in order 
to preserve the gauge invariance the field should be hermitian in Euclidean 
space-time. Then in body-fixed frame the Dirac operator of the rotating soliton 
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D{U) can be expressed in the form: 
D{U) = D{U) R=dr + h{U) + iA^ R^QR - i-f4lkAk R^QR . (34) 



Here, the rotation velocity matrix Vt in the Euchdean space-time is given by: 

= -r" = -z i?^ i? , (35) 
2 

where the derivative R is with respect to the Euchdean time r and as usual, 
summation over repeated indices is assumed. Similar to ref. the angular 
velocity VL is hermitian in Euclidean space-time. The collective variable VL 
will be quantized according to the canonical quantization rule (pTj) which in 
Euclidean space-time is: 

(36) 

where J" is the spin operator. 

Now we can write down the hermitian conjugate to (|3^) : 

D\U) = -dr + hitJ) -in- iAi R^QR - t-f47kAk R^QR , (37) 



As we can see from (p^, (|3y|), the operator ( p4D is non- hermitian. Hence, in 
general, the Euclidean action (|l^ will have the real and imaginary part. To 
make clear distinction we write down both parts explicitly: 

^eff ~ ^^eff + ^^eff ) (38) 

where 

^S^ff^—N^TrhgiD^D), (39) 

end 



QS^ff^'-N^TrlogiD/D^), (40) 



The imaginary part is finite and does not need regularization. In addition, any 
regularization of this part of the action would break several important features 
of the model. On the other hand, the real part is infinite and a regularization 
is necessary. We introduce the Schwinger proper-time regularization of the 
action. In general, for an operator AB~^ its proper-time regularized form 
reads [ITi 



Trlog(ifi~i) ^ - Tr / — e""^ - e""^ . (41) 
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Applying ( plD to the real part of the action (|39D we obtain: 



^S^^^ = ^ Tr j^e-"^'^ - e-"^i^"j , (42) 

The cut-off A is fixed in the meson sector to reproduce the physical pion 
decay constant f.^. We also subtract the vacuum contribution where the cor- 
responding operator: 

D, = dr + K (43) 
includes the unperturbed (vacuum) hamiltonian Hq: 

a ■ V 

K = — + PM + mop . (44) 



Thus the effective action includes the regularized real part (|39|) and finite 
imaginary part (^01). In fact, for the contribution from the imaginary part 
(|40| ) one can use the non- regularized expression ([21|). 

Now, our task is to evaluate the contribution (|3lD coming from the regular- 
ized real part (|39D. In this case the difficulty lies in that the expression (|i2| ) 
contains rather complicated operator exponent exp{—uD^D) in the inte- 
grand. However, because afterwards we take derivative over the fields and 
set Afj^ = 0, the only non-zero contribution will come from the part of the 
integrand linear in A^. We shall treat the angular velocity ~ as a small 
perturbation that is consistent with the large Nc limit philosophy behind the 
NJL model. Hence, for the electromagnetic current we can neglect terms 
and higher. However, as can be seen later the term ~ must be taken into 
account to evaluate the moment of inertia (see also ref. [0). 

In the next step we expand the integrand in (|3TD in terms of and fi. It is 
useful to separate 

D^D = dIDo + + V, (45) 

where 

dIDo = ~dl + , (46) 
and the sum of all perturbative terms we denote by V: 

V = Wo + Wi + W2 (47) 

with 

Wo = i[h, n]-i{Vt, dr}, (48) 

Wi = [h, lAiR^QR] + Ai{R^QR , Vt} - i{AiR^QR , d^} , (49) 
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and 

W2 = -t-f,jkAk{h, R^QR}+-f^-fkAk[R^QR, fi]-z747fc [AkR^QR, d^]. (50) 

Here [ , ] and { , } denote commutators and anticommutators, respectively. 
After that for the operator exponent Trexp {—uD'^D) in eq.(^) we apply an 
expansion 



A+B _ + / da e°^ B e^^"")^ 



1 1-/3 

+ J d(3 J da e"^ B e^^ B e^^"""^)^ +... (51) 


J. 

There for the operator A we take the part —uDqDq, which does not contain 
A^ and f2, and using the cyclic properties of Tr we get: 

2 1 

Tr e-"^'^ = Tr e'^^J^" -^z Tr (^e~"^.)^° + ^ J dp 





We compute the trace over the Euclidean time in the w-space, the Fourier 
conjugate to r: 

duj 



{T\f{d^,...)\T) 



27T 



f{-iuj,...) . 



(53) 



Applying the operator expansion (|5^), there are, in general, two types of terms 
in (^) - linear and bilinear in V . Obviously, in order to give a non-zero 
contribution to eq.(PTD all those terms should be linear A^. This implies that 
only (2) -terms will contribute to the first order in ([5^): 



- uTi( e-"^o^o 



6Wi 



(2) 



SA^iO, z) 



A„=0 



(54) 



One should also keep in mind that 14^1(2) includes also terms linear in VL. 



The terms contributing to the second order in (|5l|) have a more complicated 
structure: 



SA^iO) 



+ Wo 



SWi(2)_ 

6A„ 



(55) 



In the case that the integral over R includes both and Dab we are left with 
an additional integral over Euclidean time which includes their time-ordered 
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product dH): 



J UJ — UJ 

+ T^5{u;-J){J,, Dab], (56) 
where in r.h.s. V means a principal value. 

To summarize, in the scheme presented above the matrix element of the cur- 
rent includes non-regularized as well as regularized contributions. In gen- 
eral, there are two types of non-regularized contributions: the valence parts 
(for all form factors) and the sea parts coming from the imaginary part of the 
action. To compute those terms we use the non- regularized expression (pT]). 
For the contribution from the real part of the action we follow the expressions 
which include an explicit proper-time regulator (^T). 



To begin with explicit evaluation of the particular current matrix elements 
and the corresponding form factors we will sketch briefly the derivation of the 
Dirac sea contribution to the moment of inertia within the present scheme. 
The valence part which does not need regularization and can be found in 
ref. §: 

Qvai ^ ^ {val\Ta\n){n\n\val) ^^^^ 

Here \n) and e„ are the eigenfunctions and eigenenergies of the hamiltonian h 
and val stands for the valence level. 

The Dirac sea part of the moment of inertia originates from the regularized 
real part ^Sefj of the effective action (^) with electromagnetic field set to 
zero. The imaginary part does not contribute. In first order operator expansion 



(51) only the term contributes 



1 



40r ^ •/ 2 

{n\Ta\m)(m\n\n) , (58) 



where for convenience we introduce the complete set \m){m\ and symmetrize 
with respect to m. The second order term in expansion (|5lD includes 

only Wq. 
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/ duu / d/3 Trf e-"(i-'3)°J°o ^ ^-^f^^l^^ ^ + a ^ b V (59) 

J J V OiZa OiZu / 

1 Q 

It is easy to be seen that the above expression is symmetric in a ^ h . After 
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some straightforward calculations one gets: 

{m\Ta\n){n\,n\m) . (60) 



A^c ^ 7" du e " — e " 



9^1 /iT p2 _ p2 

^ n,m , V ^ 



- + - em - en 

.M 2 



Combining ( |55D and (0) the final result is: 



©r = Y ^ {n\Ta\m){m\n\n)n%{e^, e„). (61) 



Here is the regularization function: 
"^9 (em 5 e^) 



2 2 2 2 

du / I e~"^" — e~"^™ e„ e^"^" + e^'"*^'" 



4A/7r y v^Vu e^-e2 e^ + e„ 

(62) 

which coincides with ref. It is easy to see that the moment of inertia is 
diagonal: 

Qab = SabQ . (63) 



3 Computing the form factors 



The nucleon electromagnetic Sachs form factors are related to the matrix 
element of the electromagnetic current by: 

{N{p')\MO)mp))=GEiq') , (64) 
{N{p')\jm\N{p)) = ^ GM{q') ^ {N\{a x q)\N) , (65) 

where \N{p)) is the nucleon state of a four momentum p and proper spin 
and isospin quantum numbers. The four momentum transfer q = p' — p. The 
electromagnetic current ./^(x) is in the Minkowski space-time and Ai^ is the 
nucleon mass. The expressions for the matrix element of the current have been 
derived in the previous Section. 

The isoscalar and isovector parts of the form factors are defined by: 

Ge{m) = - G%Jm) + T3 G%jIj-^ . (66) 



The isoscalar form factors (electric and magnetic) come from the scalar part 
of the quark current matrix Q (~ 1) while the isovector form factors is from 
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the triplet part (~ r^). Since the imaginary part of the effective action in 
Euchdean space-time (|38|) is due to the isoscalar part of the charge matrix Q 
one should expect that for both isoscalar (electric and magnetic) form factors 
the Dirac sea contribution originates from the imaginary part of the effective 
action. Indeed, as can be checked directly the contribution from the regularized 
real part of the action to both isoscalar form factors is exactly zero. This means 
that for the calculation of these form factors the non-regularized expression 



([21| ) can be used. 



We start with the electric scalar form factor. Since only the isoscalar part of 
the operator Q contributes (Jo = |l) the rotational matrices R{t) vanish from 
the action. Hence, the DTZ integration in ( pi]) becomes trivial. It can be seen 



easily from (Q) and (^) using the explicit form of Wq and Wi that because 
of symmetry considerations there is no contribution from the real part of the 
action and the Dirac sea contribution comes only from the imaginary part (^Ol). 
We start with the non- regularized eg . (^1]) , expand the effective action and the 
quark propagator in f2, (p2|) -(p3|), and insert the spectral representation (pQl). 
In this case the contribution from the terms linear in f2 is exactly zero. Finally 
we arrive at: 



Y: sign(e„) ^U^) $„(£) } , (67) 



1 

~ 2 

The first term in eq.(|67|) comes from the valence quarks and the second term 
is due to the quark vacuum polarization. 

At = the isoscalar electric form factor reproduces the baryon number 
B = 1: it easy to see that the contribution from the sea vanishes exactly 
whereas the valence term gives one (normalization of ^vai)- 

As the next step, we compute the electric isovector form factor. Using the 
formulae (|2T| ) - (|29|) the evaluation of the valence part is straightforward (see 
ref. [H) and one gets 

(Ar(p') l^t^^i 

^ J ^^^^^ em Cval 

X (1/2, JgTg |{ J, , DaJI 1/2, JsTg) . (68) 

The anticommutator {J^ , -Ds;,} appears due to the symmetry of the r matrix 
element product under n ^ m . 

The contribution of the Dirac sea comes from the real part of the effective 
action and a regularization is necessary. Contributions to the first order in 
expansion (|52D come only from the term of Wi containing both A4, and VL: 



15 



Tr( e-"^o 



Do 



5Wi 



5^4(0, z) 



Nc ^ 7 du e-"^" + e-"^' 



A4=o 4v^„,„y 2 



r, |m) {<^l{zy<^^{z) ) (1/2, J3T3 |{ Jc , D,,}\ 1/2, J3T3) • (69) 



The other two terms in Wi, i[h , A^R^^^] —iiA^R^^R^ ^r} , give no 
contribution. 

The second order in (^) includes terms from Wq and Wi. 



^Jdp{l-p)Jdu M Tr(^e-"^'5^o e- 

6Wi 



6Ai{0,z) 



6A, 



A4=0 

(70) 



Because of the symmetry of the r matrix element product under n ^ m only 
the symmetric part of the r.h.s. of the integral (^) survives. After integrating 
over uj and (3 (the integrals are the same like in the case of the moment of 
inertia) one obtains: 



40F0 ^{ 4 - « 2 

%l{z) T, $„(?) ) (1/2, J3T3 \{l , Dst}\ 1/2, J3T3) . (71) 



Finally, using the quantization rule (^) the collective matrix element can be 
easily calculated: 

(l/2,4T^|{J,,D„,}|l/2,J3T3) = -l5e.(r%, 5,,, . (72) 

6 33 •^3'^3 

Combining (|68D - (^) we arrive at the following expression for the electric 
isovector form factor in Minkowski space-time: 



(*L,(-')^°*«(?)) {n\T'\val) 



+ E Ki^n., en) ($J,(i')r"$„(i')) (m Ir-^l n) 



(73) 



where the regularization function TZq is exactly the same as in the case of 
moment of inertia. Hence, formula ([73D for = is equivalent to the quotient 
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of the moment of inertia by itself: 

Gr(g^ = 0)=l, (74) 

as it should be. 

Now we proceed to compute the magnetic form factors. We start with the 
isoscalar one. In this case, the quark charge matrix Q contributes via the 
isoscalar term, ~ | 1, diagonal in the isospin space, the rotation matrices R{t) 
cancel each other and the current is ~ |7o7fe- Similar to the isoscalar electric 
form factor, the Dirac sea contribution from the real part of the action is zero: 
the term from the first order in (|5^) is proportional to (n |7o7fc| n) = and in 
the second order the terms are odd in uj and after integration vanish. Hence, 
the Dirac sea contribution comes solely from the imaginary part of the action 
and the current matrix element can be calculated directly from (pH). However, 
since the current does not include rotation matrices R{t) the first (zero order 
in VL) term from (|23| ) gives no contributions. In the leading order (linear in VL) 
the collective part includes only f2: 

(/V(p') |*t,„,,*| N(p)) = / d'. e-'^- Y. (^^('>°^^*"'^'')<" l^°l "') 

n>val 
m<.val 

x(l/2, J3T3IJCII/2, J3T3). (75) 

From (|65|) and ([75|) resolving the spin structure we get for the magnetic 
isoscalar form factor the following expression: 

^, /'^3_-iqz. qj V- *m(z)7o7k$n(z)(n|r^|m) 



Q n>val 
m<val 

(76) 

where, as usually, summation over repeated indices k,j,a is assumed. 

For the magnetic isovector form factor the current includes the isotriplet 
part of the quark charge matrix Q. Because of this we cannot get rid of the 
rotation matrix R{t) in the action. The valence part is obtained from eq.(^T]). 



Expanding the quark propagator in f2 and using the spectral representation 
([29| ) it is straightforward to calculate the contributions from the first as well 
as from the second order terms in expansion (p3|). For the first order (zero in 
Q) there is no problem with time-ordering and the collective integral ( pSj) is 
reduced to an ordinary integral of three Wigner D-functions: 

(1/2, J',n I D,, I 1/2, J3T3) = (r")^,^ {r")^,^ (77) 

6 -I3J3 J3J3 

In the second order (linear in Q) we have a time-ordered product (pHj). In 
contrast to the isoscalar magnetic form factor the quark matrix element 



17 



(n |7o7fcT^| m){m \ T°'\n) is asymmetric under n ^ m and we end up with a 
commutator in the collective part: 

(1/2, I [4 D,,] I 1/2, JsTg) = ze^"'^ (1/2, 4r^|D,,| 1/2, JsTa) (78) 

Summing up the first and the second order terms we get for the valence part: 



(iV(p')|*S7fey^|iV(p)).a/ = 



d^z e"'^-^ 



($L(^)7o7fcr"$.aK^))(l/2,J3T3|D36|l/2,J3T3) + ^ E sign(eO 



l/2,JsTs\[l,Dsb]\l/2,.m)\. (79) 



X (1/ 

The Dirac sea part originates from the regularized real part of the effective 
action (38) (in Euclidean space-time). In the first order in (^) the only term 
from W2 (|50|), which gives non-zero contribution, is i'~f4'yk{h, K^QR}. The 
other two, ^747^ [R^QR , dr] and — 747A; [R^QR , Q] simply vanish. Using ([77|) 
one gets 



00 

N,J du Tr(^e-"^oA, 



7^ 



...... / -NcT.<nien) 

[^i{z)^^-1ur^^n{z)) (1/2, J3r3 1^361 1/2, 73^3) 



^0) 



with 71^1 being the corresponding regulator: 



7^ 



Ml(^n) 



du 



20Fi yil^" 



[81) 



To the second order 
1 



4 



j d(3 j duu Tr(^ e-"(i-^)^J^o 6""^^^ 

5Wl 



Do 



Mfc(0,f) 



f82) 



contribute i747A: {/i, R^QR} and 2747^ [-R^Q-R, <9r] • Using the asymmetry of 
the quark matrix element (n 1 747^x^1 m)(m |r"| n) under n ^ m one can 



show that the symmetric part of the r.h.s. of the integral (|56D vanishes. Only 
the asymmetric part survives: 



/ d/? / duu / / 

29 y 7 J 271 J 271 



du f du' 



1 



UJ — UJ' 
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x($]„(z)747A:T''^>n(^))(^|r1m)(l/2,J3T3|[4D,fc]|l/2, J3T3) (83) 

Using ( |78|) and performing the integrations over uo and u we get for the second 
order contribution in the form: 

E^M2(em,en)($L(i')7o7fcrb$„(i'))(n |r^| m)(l/2, J3T3 |[J,, 1/2, J3T3) , 

(84) 

with a second regulator 7^^2- 



1 



1 







Summing all obtained contributions (|79D -(^4D and using eqs.(l65|) and (|66D, 
we obtain the final expression for the magnetic isovector form factor in the 
following form: 



2® n^val ^" 



+ 4^"'''E ^M2(en.,en) ($U^l7o7fc <f n^(i')) |r» | . (86) 



The proton and neutron form factors are defined, respectively, as sum and 
difference of the isoscalar and isovector form factors: 

r<P — 1 [ r'T=0 I ryT=l 1 
'^E,M — 2 [^E,M ~^ ^E,M i ; 

_ 1 r r~iT=0 _ /^T=l 1 
'^E,M — 2 l^E,M '^E,M\ ; 

and the numerical results for them are presented in Section 4. 




4 Numerical results 



To compute observables we use and numerical method of Ripka and Kahana 



T8| for solving the eigenvalue problem in finite quasi-discrete basis. We con- 



sider a spherical box of a large radius D and the basis is made discrete by 
imposing a boundary condition at r = Also, it is made finite by restricting 
momenta of the basis states to be smaller than the numerical cut-off K^ax- 
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Both quantities have no physical meaning and the results should not depend 
on them. The typical values used are D ~ 20/M and K^ax ~ 7M. 

In addition, all checks concerning the numerical stability of the solution with 
respect to varying box size and choice of the numerical cut-off have been done 
and the actual calculation is completely under control. 

The actual calculations are done by fixing the parameters in meson sector in 
the well known way to have /^r = 93 MeV and = 139.6 MeV . This 
leaves the constituent quark mass M as the only free parameter. 

The proton and neutron electric and magnetic form factors are displayed in 
Figs.0-|^. The theoretical curves resulting from the model are given for the 
following four values the constituent quark mass: 370,400,420 and 450 MeV. 
The magnetic form factors are normalized to the experimental values of the 
corresponding magnetic moments at = 0. With one exception of the neutron 
electric form factor (Fig.^), all other form factors fit to the experimental data 
quite well. The best fit is for the constituent quark mass around 400-420 MeV. 



Figs.|Hl| 



As can be seen the only form factor which deviates from the experimental 
data is the neutron electric form factor and this requires some explanation. 
Obviously, this form factor is the most sensitive for numerical errors. According 
to the formula ( p8D the form factor has been calculated as a difference of 
the electric isoscalar and electric isovector form factors that were the direct 
output of our code. Both form factors were of order 1.00 and calculated by 
the code with a high enough accuracy. However, the resulting neutron form 
factor has the proper (experimental) value of order 0.04, i.e. about 4% of the 
value of its components. This means that a small numerical error for one of 
the components can be enhanced by a factor 50. Hence, small numerical errors 
together with the applied approximations (like the l/N^ approximation behind 
the model and neglecting the boson-loop effects as well as higher order fermion 
loops) are strongly magnified resulting in a deviation from the experimental 
data for momentum transfers above 100 MeV/c. 



Table |l] 



As the next step, one can compute other electromagnetic observables: the 
mean squared radii, the magnetic moments and the nucleon-A splitting. The 
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static nucleoli properties, in particular the charge radii and the magnetic mo- 
ments can be obtained from the form factors: 



^T^0,1 ^ ^T^0,1(^,^ ^^^^ _ (go) 

For the quark masses 370, 420 and 450 MeV and pion mass m^r = 14QMeV 
the calculated values are presented in Table |l|. The values for the axial vector 



(89) 

g2=0 



coupling constant qa-, presented for completeness, are from ref. [|12|. For com- 
parison, in Table |^ we give the theoretical values of the same quantities but 
with the physical pion mass set to zero. 

The chiral limit (m^ 0) mostly influences the isovector charge radius. In- 
deed, as it should be expected |[19|,|6[|, in our calculations with zero pion mass 



the isovector charge radius diverges in chiral limit. It is illustrated in Fig.|^, 
where the isovector charge radius is plotted vs. the box size D. The Dirac sea 
contribution to the radius grows linearly with D and diverges as -D oo. 
Because of this that quantity (and the relative quantities) is not included in 
Table H. The other observable strongly influenced by the chiral limit is the 
neutron electric form factor. For the m,r ^ the discrepancy from the exper- 
iment is by almost a factor two larger than in the case ^ 0. The other 
observables differ in the chiral limit by about 30%. The comparison of the 
values in the two tables indicate that taking the physical pion mass gives us 
a best fit with a much better agreement with the experimental data. To be 
particular, in chiral limit (ttItt = the calculations for the electric observables 
suggest a high value (M ~ 450 MeV) while for the magnetic ones indicate 
M ~ 370 MeV which is not the case with physical pion mass. In addition, in 
chiral limit we observe much larger contribution from the sea effects, about 
50% of the total value. 



Table g 



The results of Table (m^ = 140 MeV) again indicate the value ~ 400—420 
MeV for the constituent quark mass, in agreement with the conclusion drawn 



from the form factor curves. The same has been suggested earlier P,|TT 
where smaller number of observables has been taken into account. With the 
exception of the neutron electric squared radii, to which remarks similar to 
the case of the neutron electric form factor are applicable, the contribution of 
the valence quarks is dominant. However, the contribution of the Dirac sea is 
non- negligible, it lies within the range 15 - 40%. 
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One can notice, that the numerical results for the nucleon-A mass splitting 
(Ma — Mtv), the mean squared proton, isoscalar and isovector electric radii 
and the axial coupling constant (gA), as well as the proton electric and mag- 
netic and neutron magnetic form factor, differ from the experimental data by 
no more then about ±5%. Finally, for the magnetic moments we have got 
results 20-25% below their experimental values. Despite of this underestima- 
tion of both magnetic moments we have found surprisingly good result for 
the ratio fip/ fin which is far better than in other models. The agreement with 
the experimental value is better than 3% for the constituent quark mass 420 
MeV. So good result has been obtained neither in the Skyrme model nor in 
the linear sigma model. 

In section 3 we have found that the magnetic isovector form factor is the only 
one which includes non-zero contributions in both leading (A*"^) and next to 
leading order (1/A^c) in ^- However, the numerical calculations show that the 
{1/Nc) rotational corrections do not affect the g^-dependence (slope) of the 
form factor but rather the value at the origin (G^^^(g^ = 0) = fi^^^) which 
is the isovector magnetic moment. It is not surprising since (as can be seen 
from eq. (pB]) ) in both leading and next to leading order terms the shape of the 
wave functions determines the g^-dependence of the form factor whereas 
the value at = depends on the particular matrix elements included. In 
leading order the isovector magnetic moment is strongly underestimated 0. 
As has been shown in ref. []T2[ the enhancement for this quantity due to the 
1/Nc corrections is of order {Nc + 2)/Nc and it improves considerably the 
agreement with experiment. However, as can be seen from Table |l| it is still 
below the experimental value by 25 %. 

The isoscalar and isovector electric mean squared radii are shown in Figs]^,| 
as functions of the constituent quark mass. The same plot but for the experi- 
mentally measured quantities: the proton and neutron electric charge radii, is 
given in Fig.^ In these plots the valence and sea contributions are explicitly 
given (dashed and dash-dotted lines, respectively). As could be expected, for 
the isoscalar electric charge radius the valence part is dominant (about 85%), 
due to the fact that there is no 1/Nc rotational contributions. This is not 
true for the isovector electric charge radius, where the sea part contributes to 
about 45% of the total value (Fig.|^). Also, this effect can be seen from the 
proton and neutron charge radii (Fig.^ which are linear combinations of the 
isoscalar and isovector ones. For proton the sea contribution is about 30%. 
However, for the neutron charge radius the negative sea part is dominating 
and the valence contribution is negligible. 
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Fig.^ and g 



In addition, in Fig.|TO| we plot the proton and neutron charge distribution for 
the constituent quark mass M = 420 MeV. For the proton we have a positive 
definite charge distribution completely dominated by the valence contribution 
whereas in the case of the neutron the Dirac sea is dominant. In accordance 
with the well accepted phenomenological picture one realizes a positive core 
coming from the valence quarks and a long negative tail due to the polarization 
of the Dirac sea. Using the gradient expansion the latter can be expressed in 
terms of the dynamical pion field - pion cloud. 



For completeness, in Fig.[TT| we present also the magnetic moment density 
distribution for proton and neutron for the constituent quark mass M = 420 
MeV. The sea contribution becomes non-negligible for distances greater than 
0.5 fm. Due to the relatively large tail contribution to the magnetic moments 
the sea contribution to these quantities is about 30%. 



5. Summary and conclusions 



Our numerical results support the view that the chiral quark soliton model 
of the Nambu-Jona-Lasinio type offers relatively simple but quite success- 
ful description of some low-energy QCD phenomena and, in particular, of 
the electromagnetic properties of nucleons. In parameter free calculations we 
have obtained overall good results for the electromagnetic form factors, the 
mean squared radii, the magnetic moments and the nucleon-A splitting. Us- 
ing = 93 MeV and m^r = 139.6 MeV and a constituent quark mass 
of M = 420 MeV the isoscalar and isovector electric radii are reproduced 
within 15%. The magnetic moments are underestimated by about 25%, how- 
ever, their ratio /ip//i„ is almost perfectly reproduced. The g-dependence of 
G%{q^), G\j{q^) and G^ji^q^) is very well reproduced up to momentum transfer 
of 1 GeV. The neutron form factor G^{q'^) is by a factor of two too large for 
> 100 MeV"^. One should note, however, that G%{q^) is more than an order 
of magnitude smaller than G^^{q^) and as such it is extremely sensitive to both 
the model approximations and the numerical techniques used. Here, it turns 
out that the agreement is noticeably worse if the chiral limit (m,r 0) is used. 
This can be easily understood since the pion mass determines the asymptotic 
behaviour of the pion field. Altogether we conclude that the chiral quark soli- 
ton model based on a bosonized NJL type lagrangean is quite appropriate for 
the evaluation of nucleonic electromagnetic properties. 
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Fig. 1. Diagrams corresponding to the expansion in Q of the electromagnetic current 
matrix element: a) the valence contribution and b) the Dirac sea contribution. 

Fig. 2. The proton electric form factor for the momentum transfers below 1 GeV. 

Fig. 3. The neutron electric form factor for the momentum transfers below 1 GeV. 

Fig. 4. The proton magnetic form factor normalized to the experimental value of 
the proton magnetic moment at = for the momentum transfers below 1 GeV. 
The normalization factor can be extracted from Table 1. 

Fig. 5. The neutron magnetic form factor normalized to the experimental value of 
the neutron magnetic moment at = for the momentum transfers below 1 GeV. 
The normalization factor can be extracted from Table 1. 

Fig. 6. The isovector charge radius as a function of the box size D for m.,^ = and 
m-j, = 139.6 MeV. 

Fig. 7. The isoscalar electric charge radius as a function of the constituent quark 
mass M. The valence and sea parts are marked by the dashed and dashed-dotted 
lines. 

Fig. 8. The isovector electric charge radius as a function of the constituent quark 
mass M. The valence and sea parts are marked by the dashed and dashed-dotted 
lines. 

Fig. 9. The electric charge radii of proton and neutron as functions of the constituent 
quark mass M. The valence and sea parts are marked by the dashed lines. 

Fig. 10. The charge density distribution of the proton (lower) and neutron (upper) 
for the constituent quark mass M = 420 MeV. 

Fig. 11. The magnetic moment density of proton and neutron for the constituent 
quark mass M = 420 MeV. 
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Table 1 



Nucleon observables calculated with the physical pion mass. 



Quantity 


Co 

370 MeV 


Qstituent Quark IVl 

420 MeV 


ass 

450 MeV 


Exper. 


total 


sea 


total 


sea 


total 


sea 


< >T=o ifm'^] 


0.63 


0.05 


0.52 


0.07 


0.48 


0.09 


0.62 




1.07 


0.33 


0.89 


0.41 


0.84 


0.45 


0.86 




0.85 


0.19 


0.70 


0.24 


0.66 


0.27 


0.74 


< r2 >„ [fm^] 


-0.22 


-0.14 


-0.18 


-0.17 


-0.18 


-0.18 


-0.12 


HT=o [n.m] 


0.68 


0.09 


0.62 


0.03 


0.59 


0.05 


0.88 


Ijlt=i [n.m] 


3.56 


0.77 


3.44 


0.97 


3.16 


0.80 


4.71 


fip [ri.rn.] 


2.12 


0. 13 


2.03 


0.50 


1.86 


0.13 


2.79 


fin [n.m.] 


-1.44 


-0.34 


-1.41 


-0.47 


-1.29 


-0.38 


-1.91 


\l^p/l^n\ 


1.47 




1.44 




1.44 




1.46 


Ma- MN[MeV] 


213 




280 




314 




294 


9A 


1.26 


0.08 


1.21 


0.11 


1.13 


0.06 


1.26 



Table 2 

Nucleon observables calculated with the zero pion mass. 



Quantity 


Co 

370 MeV 


nstituent Quark IVI 

420 MeV 


ass 

450 MeV 


Exper. 


total 


sea 


total 


sea 


total 


sea 


< >T=o [fm^] 


0.88 


0.20 


0.66 


0.26 


0.61 


0.23 


0.62 


fJ'T=o [n.m.] 


0.66 


0.07 


0.59 


0.09 


0.57 


0.09 


0.88 


/xt=i [n.m.] 


4.61 


1.59 


4.38 


1.89 


3.91 


1.48 


4.71 


fj,p [n.m.] 


2.63 


0.83 


2.49 


0.99 


2.24 


0.79 


2.79 


fin [n.m] 


-1.97 


-0.76 


-1.89 


-0.90 


-1.67 


-0.70 


-1.91 


\\^vli^n\ 


1.34 




1.34 




1.34 




1.46 


Ma - MN[MeV] 


221 




261 




301 




294 
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